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We studied local electronic states near Zn in op- 
timally doped YBa2(Cui_ a; Zn I )307_i and underdoped 
YBa2(Cui_ a ;Zn a; )408 via satellite signals of plane-site Cu(2) 
nuclear quadrupole resonance (NQR) spectra. From the rel- 
ative intensity of Cu NQR spectra, the satellite signals are 
assigned to Zn-neighbor Cu NQR lines. The Cu nuclear spin- 
lattice relaxation time of the satellite signal is shorter than 
that of the main signal, which indicates that the magnetic cor- 
relation is locally enhanced near Zn both for the underdoped 
and the optimally doped systems. The pure YBa2Cu40g 
is a stoichiometric, homogenous, underdoped electronic sys- 
tem; nevertheless, the Zn-induced inhomogeneous magnetic 
response in the CuC>2 plane is more marked than that of the 
optimally doped YBa2Cu307-«. 
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I. INTRODUCTION 

The effect of nonmagnetic impurity Zn in high- 
temperature cuprate superconductors has attracted great 
interests. The impurity Zn induces Curie magnetism, a 
local staggered susceptibility in the normal state [1-5], 
an in-gap antiferromagetic dynamical spin susceptibil- 
ity, and a virtual bound electronic state around Zn 
in the superconducting state [6-8]. The perturba- 
tion theory for the mean-field states close to two- 
dimensional antiferromagnetic instability allows coex- 
istence of the staggered-type, generalized Ruderman- 
Kittel-Kasuya-Yosida (RKKY) oscillation and the local 
density of states of electrons with a zero bias resonance 
peak at finite temperatures [9-11]. However, recent nu- 
merical calculation using Giitzwiller approximation indi- 
cates that the frozen staggered moments are locally re- 
duced by Zn in the superconducing state of a moderately 
doped two-dimensional t — J model [12]. According to 
this theory, the zero bias resonance peak in the local den- 
sity of states and the local staggered moments are exclu- 
sive at T=0 K. It remains to be solved whether the Zn 
strengthens [13,14] or weakens [15] local magnetic corre- 
lation around Zn in the optimally doped or overdoped 
system. 

Nuclear magnetic resonance (NMR) and nuclear 
quadrupole resonance (NQR) are powerful techniques 
to study the local states on crystalline imperfec- 
tion, dilute exchange spin networks, Friedel oscilla- 



tion, RKKY oscillation, or Kondo screening effect 
[16-20]. We focus on the local electronic states 
near Zn in YBa 2 (Cui_ a: Zn x )30 t/ (y =7-5) (Y123) and 
YBa 2 (Cui_ a; Zn 2: )40 8 (Y124) via Zn-induced Cu NQR 
satellite spectra [13,21]. Y123 is a nearly optimally doped 
system (oxygen content y ~6.9), while Y124 is a stoichio- 
metric underdoped one. Zn-neighbor 89 Y NMR signals 
and their relaxation times have intensively been stud- 
ied for Y123 [3] and for Y124 [22]. However, there 
is only a few direct observation of the local electronic 
state via Zn-neighbor Cu NQR/NMR signals; a report 
on the satellite Cu nuclear spin-lattice relaxation of a 
non-superconducting Y124 with heavily Zn doping [21] 
but no report on the superconducting Y123 or Y124 with 
small concentrations of Zn. 

In this paper, we report on a study of satellite signals 
in the planar Cu(2) NQR spectra of Zn-substituted su- 
perconductors Y123 and Y124. The wipeout effect on 
the Cu NQR main signals and the growth of the Zn- 
induced Cu satellite signals indicate that the satellite but 
not the main signal includes the Zn-neighbor Cu nuclei. 
We demonstrate that the satellite Cu nuclear spin-lattice 
relaxation time is shorter than the main one in the super- 
conducting Y123 and Y124. The local antiferromagnetic 
correlation near Zn is enhanced even for the optimally 
doped system. 

II. EXPERIMENTAL 

High quality powder samples of a pure Y123 (y=6.98 
studied in Ref. [23]), a slightly oxygen-deficient Y123 
(y=6.89 in Refs. [24,13]), Zn-substituted Y123 (the an- 
nealed samples in [25]), a pure Y124 and Zn-substituted 
Y124 [8,26] were prepared for the present study. The su- 
perconducting transition temperatures are T c ~92 K both 
for Y123 with y=6.98 and with y=6.89; T c =81, 69, and 
46 K for Y123 with Zn content x=0.007, 0.017, 0.033; 
T c =82 K for the pure Y124, and T c =68, 56, 37, and 27 K 
for Y124 with Zn content £=0.005, 0.01, 0.015, 0.022, re- 
spectively. To estimate the relative intensity of Cu NQR 
spectra, the respective samples of Y123 and of Y124 with 
nearly the same volume were prepared. The oxygen con- 
tents conform to the values in the previous studies. 

A phase-coherent-type pulsed spectrometer was uti- 
lized to perform Cu NQR experiments. The zero-field Cu 
NQR frequency spectra with quadrature detection were 
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obtained by integration of Cu nuclear spin-echoes with 
changing the frequency in a two-pulse sequence (tt/2- 
t-7t echo). A typical width of the first exciting n/2- 
pulse t w was about 3 /is (the excited frequency region 
v\ ^83 kHz from 2-KV\t w =-K /2). The Cu nuclear spin- 
lattice relaxation curves were measured by an inversion 
recovery technique, where the nuclear spin-echo M(t) was 
monitored as a function of the time t in a sequence of 
7T — t — 7r/2 — t — 7T echo. The pulse interval for the NQR 
spectrum measurement was t=15 fis. 



III. RESULTS AND DISCUSSION 
A. Main and satellite 63 65 Cu(2) NQR spectrum 

1. Frequency- distributed Cu(2) nuclear spin-lattice relaxation 

Figure 1 shows 63,65 Cu(2) NQR spectra (upper pan- 
els) and the 63 Cu(2) nuclear spin-lattice relaxation curves 
p(t) = 1 — M(t)/M (oo) (lower panels) for the oxygen- 
deficient Y123 (left), Zn-substituted, nearly optimally 
doped Y123 (middle), and Zn-substituted, underdoped 
Y124 (right) at T=4.2 K. The Cu NQR spectra are nor- 
malized by the respective values of peak intensity. Zn 
doping induces lower frequency tails or satellite signals 
both for Y123 and Y124, which is known in Refs. [13,21]. 
Oxygen deficiency is also known to yield similar low fre- 
quency satellite or shoulder in the Cu(2) NQR spectrum 
for Y123 without Zn [27-30]. Since Y124 has no appre- 
ciable oxygen deficiency, the satellite signals result from 
purely Zn doping effect. 

For oxygen deficient Y123 in the lower panel of Fig.l, 
the satellite and the main signals show nearly the same 
relaxation curve p(t), and their relaxation times T\ are 
longer than a few seconds. This is consistent with the 
previous reports [27,30]. On the other hand, for Zn- 
doped Y123 and Y124, the Zn-induced satellite signal 
shows a shorter relaxation time than the main signal, 
and their relaxation times are shorter than 1.0 s. Thus, 
the lower frequency satellite in the Cu NQR spectrum 
for Zn-doped Y123 possesses an electronic state different 
from that for the slightly oxygen deficient one. 

Our observation for Zn-substituted Y123 and Y124 
that the lower frequency signals in the broad NQR spec- 
tra show the shorter T\ is in parallel to that for Zn-free 
underdoped La 2 -xSr x Cu0 4 [31,32]. For Y124 (Y123), 
we confirm that the satellite recovery curve p(t) is nearly 
independent of frequency over 28.4-28.6 MHz (30.1-30.5 
MHz) and so is the main one over 29.2-30.2 MHz (31.5- 
31.6 MHz). Here, the excited frequency region by the 
pulse is v\ ~83 kHz at each NQR frequency. Thus, 
the frequency dependence of the relaxation curve is not 
monotonic. 
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FIG. 1. Upper panels: Zero-field frequency spectra of 
the planar 63 ' 65 Cu(2) NQR for oxygen-deficient Y123 with 
y=6.89 and j/=6.98 (left), Zn-substituted Y123 (middle), and 
Zn-substituted Y124 (right) at T=4.2 K. Lower panels: The 
main and the satellite 63 Cu(2) nuclear spin-lattice relaxation 
curves p(t) for the representative samples at T=4.2 K. The 
solid curves are fitted by theoretical recovery curves (see 
the text). For Y124 (Y123), the satellite recovery curve 
p(t) is nearly independent of frequency over 28.4-28.6 MHz 
(30.1-30.5 MHz) and so is the main one over 29.2-30.2 MHz 
(31.5-31.6 MHz). 



2. Isotope dependence of Cu(2) nuclear spin-lattice 
relaxation 

Figure 2(a) shows the 63 < 65 Cu(2) NQR spectrum of 
Y123 with Zn content of £=0.033 (T c =46 K) at T=4.2 
K and the hot spot regions (shaded bars) excited in this 
spectrum to measure the isotope dependence of Cu(2) 
nuclear spin-lattice relaxation curves. Figure 2(b) shows 
the respective 63 ' 65 Cu(2) nuclear spin-lattice relaxation 
curves p(t) for the above Y123 at T=4.2 K. The dashed 
and the solid curves are the respective least-squares fit- 
ting results for 63 Cu(2) and 65 Cu(2), using a theoretical 
function including a stretched exponential function (see 
below). The nuclear spin-lattice relaxation rate due to 
a magnetic origin is proportional to the square of the 
63 ' 65 Cu nuclear gyromagnetic ratio ( 63 ' 65 7„) 2 . Since we 
could not find any large deviation from the scaling with 
respect to the reduced time ( 63 ' 65 7 n ) 2 <, the overlapping 
effect of the broad 65 < 63 Cu(2) lines is small. 



B. Zn-substitution effect on 63 Cu(2) NQR spectra 

1. Zn-induced wipeout effect and satellite 

Figure 3 shows the 63 Cu(2) NQR spectra with Zn dop- 
ing at T=100 K(> T c ) for Y123 with Zn contents of x=0, 
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FIG. 2. (a) The 63 ' 65 Cu(2) NQR spectrum of Y123 with 
Zn £=0.033 (T c =46 K) at T=4.2 K and the hot spot regions 
(shaded bars) excited by rf pulses to measure the nuclear 
spin-lattice relaxation, (b) The 63 ' 65 Cu(2) nuclear spin- lattice 
relaxation curves p(t) measured at the respective frequencies 
denoted in (a). The dashed ( 63 Cu(2)) and the solid ( 65 Cu(2)) 
curves are the least-squares fitting results using a theoretical 
function including a stretched exponential function (see the 
text). 



0.007, 0.017, 0.033 (a), and for Y124 with Zn x=0, 0.01, 
0.015, 0.022 (b) from the top to the bottom. The normal- 
state Cu NQR intensity is not affected by the Meissner 
shielding effect. The main signal intensity of the 63 Cu(2) 
NQR spectrum around 31.5 MHz for Y123 and around 
29.8 MHz for Y124 rapidly decreases with increasing Zn 
concentration. Alternatively, the lower frequency tail or 
the broad satellite signal increases with increasing Zn 
concentration. The solid and the dotted curves are nu- 
merical simulations using three Gaussian functions (Mi, 
M 2 , and S). By using two Gaussian functions, Mi and 
M2, we can reproduce the main spectra. The shaded 
Gaussian S is the satellite spectrum. 



2. Relative intensity 

In Figs. 3(c) and 3(d), the estimated intensity of the 
main spectrum / ma in(=Mi+M 2 ) (upward closed trian- 
gles) and that of the satellite spectrum I sate (=S) (open 
circles) normalized by 1(0) are plotted as functions of the 
plane-site Zn(2) concentrations c=3x/2 for Y123 (c) and 
of c=2x for Y124 (d). The open squares are the total 
intensity (=M1+M2+S). The values of I m ain an d hate 
are estimated after T 2 correction using M(0) from fitting 
M(r > 15^s)=M(0)cxp[-2r/r 2L - 0.5(2t/T 2 g) 2 ] with 
the fitting parameters M(0), T 2£ and T 2 g- 

The dashed curves are numerical calculations of 
Cu intensity with various Zn configurations, using 
the binomial-probability n ^Bk(c)= n ^CkC k (l — c) n ^~ k , 
where n(j) is the number of the jth nearest neighbor 
(nn) Cu sites. (1-c) is the probability that a plane site 
is occupied by a Cu atom but not a Zn atom, or it 
is the total Cu NQR intensity with Zn doping. With 
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FIG. 3. Zn-doping effect on the plane-site 63 Cu(2) NQR 
spectra for Y123 (a) and for Y124 (b) at r=15 /is and at 
T=100 K (> T c ). The line shapes are scaled before T2 cor- 
rections. The relative intensity Imain/ 1(0) (upward closed 
triangles), I sa te/I(0) (open circles), and the total intensity 
(open squares) estimated after T2 corrections are plotted as 
functions of the in-plane Zn(2) concentration c(=3a;/2) for 
Y123 (c) and (=2x) for Y124 (d). The dashed curves in (c) 
and (d) are 7(c)/7(0)=(l - c) 5 (l), (1 - c) 13 (2), (1 - c) 21 (3), 
8c(l - c) 8 (4), and 4c(l - c) 4 (5). 
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FIG. 4. Schematic illustrations of (a) a 63 Cu(2) NQR spec- 
trum and (b) a top view of the Cu02 plane with a Zn impurity. 
The shaded area is assigned to the Zn-neighbor Cu sites. 



nn Cu sites to Zn are unobservablc. The Zn-nearest- 
neighbor Cu NQR spectrum will be shifted and broad- 
ened due to local crystalline or electric deformation, be- 
ing unobservable probably due to a shorter relaxation 
time than a few fis. For Y123, the satellite signal agrees 
with either of the 4th nn Cu site or the 2nd and the 
3rd nn Cu site to Zn [the case (4)]. With x <0.017, the 
former assignment is consistent with those to the total 
and the main signals, whereas with x=0.033, the latter 
is consistent with those to the total and the main ones. 
For Y124, the satellite signal agrees with the 1st, 2nd, 
or 3rd nn Cu site to Zn [the case (5)]. However, this is 
inconsistent with the up-to 3rd nn wipcout effect on the 
total intensity. Thus, the satellite resonance of Y124 may 
arise from the fraction of the 4th nn Cu site to Zn. 



respect to Zn configuration, the probability of finding 
Cu atoms can be seen in the decomposition of (l-c)=(l- 
c )Ilj=iJ2k=on(j)Bk( c )- In Fig. 4, schematic illustrations 
of a 63 Cu(2) NQR spectrum (a) and the Cu02 plane with 
a Zn impurity (b) are given for a guide to the calculation. 

The calculated curves on the Cu NQR intensity are 
in what follows. In Figs. 3(c) and 3(d), the decreas- 
ing functions with increasing c are 1(c)/ 1 (0)= (1 — 
c)Mc)=(l - cf (1), (1 - c)( 4 i? ( C )) 3 = (l - c) 13 (2), 
and (1 - c)( 4 B (c)) 4 =(l - c) 21 (3). These are the proba- 
bility of finding the Cu atoms with all Cu atoms in their 
1st nn positions (1), in their 12 neighboring (up-to 3rd 
nn) positions (2), and in their 20 neighboring (up-to 4th 
nn) positions (3). In other words, the diminished Cu 
NQR signals come from the Cu atoms with at least one 
Zn atom in their 1st nn positions (the 1st nn wipcout 
effect) (1), in their up-to 3rd nn positions (up-to 3rd nn 
wipeout effect) (2), and in their up-to 4th nn positions 
(up-to 4th nn wipeout effect) (3). In Figs. 3(c) and 3(d), 
the increasing functions with increasing c are I (c) / 1 (0)= 
(l-c) 8 B 1 (c)=8c(l-c) 8 (4), and (1 -c) 4 S 1 (c)=4c(l -cf 
(5). These are the probability of finding the Cu atoms 
with one Zn atom in their 4th nn positions or in the 2nd 
and the 3rd nn positions (4), and in their j(=l, 2, 3)th 
nn positions (5). 

The relative intensity of the satellite Cu NQR spec- 
trum of Y124 in Fig. 3(d) is smaller than that of Y123 
in Fig. 3(c). For Zn-substituted L^-zSr^CuO^, the dis- 
appearance of the Cu NQR signal without any satellite 
signal is observed [33]. Thus, it is likely that the observ- 
able satellite signal depends on the carrier-doping level 
or the underlying magnetic correlation. 

As seen in Figs. 3(c) and 3(d), the Zn-doping depen- 
dence of the total intensity for Y123 (x < 0.0 17) and Y124 
(x <0.022) agrees with that of the case (3) within the ex- 
perimental accuracy, where the up-to 3rd nn Cu sites to 
Zn are unobservable. The Zn doping dependence of the 
main intensity for Y123 (x <0.017) and Y124 (x <0.022) 
agrees with that of the case (2), where the up-to 4th 



3. "All-or-nothing" model 

In the above analysis, one should note that the 
so-called "all-or-nothing" model is implicitly assumed, 
where some nuclei are fully observable but the others are 
fully unobservable [37]. This enables us to estimate the 
number of observed nuclei from the observed NQR signal 
intensity. In the actual sample, however, there must be 
T2 distribution. Figure 5 illustrates a two-pulse sequence 
with a Cu NQR spin-echo in a rotating frame with a Cu 
NQR frequency in the "all-or-nothing" model (a) and in 
an actual system with T2 distribution (b). In the case 
of Fig. 5(b), some of the observable nuclear spin-echoes 
after the two pulses will be underestimated due to a some- 
what shorter T2 than the totally corrected T2 decay curve 
exp(-2T/T2), even when the Gaussian-times-exponential 
function empirically well describes M(r). Conversely, the 
I z fluctuation effect due to a short T\ [34] or the motional 
narrowing effect [35,36] on the transverse nuclear spin re- 
laxation causes an exponential decay in a longer delay 
time r but in a Gaussian decay in a shorter time r, then 
the total correction using a simple exponential function 
within the observable r range will lead to overestimation 
of the extrapolated M(0). For both the main and the 
satellite signals, such a T 2 distribution or dynamical ef- 
fect can cause overestimation or underestimation of the 
intensity. Thus, one should be careful with the results 
from intensity analysis. One can safely mention that the 
main signal does not come from the Cu sites near Zn but 
that the satellite one from the Zn-neighbor Cu sites. 

4- Cu NQR frequency 

Here, we briefly mention the Cu NQR frequency The 
Cu NQR frequencies of stoichiometic cuprate oxides have 
been analyzed with several models [38]. For non- 
stoichiometic or impure compounds, although there is 
an attempt with numerical simulation [32], it is more 
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FIG. 5. Illustrations of a pulse sequence of the first ex- 
citing 7r/2-pulse, the refocusing pulse 7r-pulse, and an NQR 
spin-echo in a rotating frame with an NQR frequency in the 
"all-or-nothing" model (a) and in an actual system with T2 
distribution (b). The "40%-visible" nuclei are supposed to 
have a shorter T2 component with 60 % than the totally cor- 
rected T2 decay curve. The solid and the dashed decay curves 
are transverse relaxation functions. 
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FIG. 6. Zn-doping effect on the main 63 Cu nuclear 
spin-echo recovery curves p(t) for Y123 with :r=0.017 (a) 
and 0.033 (b) at T=100 K. The solid curves are the 
least-squares fitting results using a theoretical function in- 
cluding a stretched exponential function (see the text). 



complicated, because local crystal strain [39] or elec- 
tric charge oscillation [40] can affect the line profile of 
Cu NQR spectrum. Thus, we do not employ a specific 
model to analyze how the satellite Cu NQR frequency 
is shifted to the lower side than the main one. We just 
emphasize that Cu 2+ is a Jahn- Teller ion but Zn 2+ is 
a non- Jahn- Teller ion, so that the local crystalline dis- 
tortion [41] must modify the electric field gradient ten- 
sor at Zn-ncighbor Cu sites, differently from that at Cu 
sites away from Zn or from that in the pure Y123 [42] or 
Y124 [43]. 
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FIG. 7. Temperature dependence of the 63 Cu nuclear 
spin-echo recovery curves p(t) at the main and the satellite 
Cu NQR frequencies for Y123 (upper) and Y124 (lower). The 
solid curves are the least-squares fitting results using a the- 
oretical function including a stretched exponential function 
(see the text). 



C. Main and satellite 63 Cu(2) nuclear spin-lattice 
relaxation 

1. Zn- substitution effect on 63 Cu(2) nuclear spin-lattice 
relaxation curves 

Figure 6 shows the main 63 Cu(2) nuclear spin- 
lattice relaxation curve for Y123 with Zn concentrations 
x=0.017 (a) and 0.033 (b) at T=100 K. The result for 
Y124 has been reported in Rcf. [26]. For Y123 with a 
small Zn concentration of x=0.017, the recovery curve 
p(t) relaxes more quickly than that for pure Y123, which 
is consistent with the result for Y124 with Zn x <0.010. 
With further doping of x=0.033, p(t > 1.0 ms) relaxes 
more slowly than that for pure Y123, which is similar to 
the result for Y124 with heavy doping of Zn x=0.022 at 
T=200 K [26]. 

Figure 7 shows temperature dependence of the main 
and the satellite 63 Cu nuclear spin-lattice relaxation 
curves p(t) for Zn-substituted Y123 with x=0.017 (up- 
per) and Y124 with x=0.01 (lower). For Y123, the satel- 
lite signal recovers more quickly than the main one below 
T c , but they show nearly the same recovery curves above 
T c . For Y124, however, both above and below T c , the 
satellite signal recovers more quickly than the main one. 
With further Zn doping for Y124 (x=0.022), the differ- 
ence between the main and the satellite signals tends to 
be reduced (not presented here), being consistent with 
the non-superconducting Y124 [21]. The result in Fig. 
7 is evidence for different magnetic correlation between 
the main and the satellite signals in the superconductors 
with small concentrations of Zn. 
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2. An analysis of nonexponential relaxation curves 

In Figs. 1, 2, 6, and 7, all the recovery curves are non- 
exponential functions. If the satellite signal comes from 
the same j th nn Cu sites by Zn and if the position depen- 
dence of Zn-induced relaxation process, i.e. the spatial 
distribution of Zn-induced moments or spin density, is in 
axial symmetry along the c axis, then the recovery curve 
must be a single exponential function. However, the ac- 
tual curves are nonexponential ones. Thus, the satellite 
signals consist of non-single Cu sites, or the Zn-induced 
moment or spin density shows non-axial symmetry, e.g. 
d-wave symmetry. 

Nonexponential relaxation is frequently observed in 
the impure materials. We have applied a conventional 
impurity-induced NMR relaxation theory [44] to the im- 
pure high-T c cuprate superconductors and quantum spin 
systems [8,26,45-47], which has been developed for di- 
lute Hciscnbcrg spin systems and dilute magnetic: al- 
loys. We believe that the analysis yields a vital, phys- 
ical framework on the strongly correlated impure sys- 
tems. The experimental p(t) is analyzed by the expo- 
nential function times a stretched exponential function 
p(t) = p(0) exp[-M/(Ti)ho* ~ y/3t/n]. p(0), (Ti) hort 
and T\ are the fitting parameters. The numerical factor 
3 is introduced to conform to the conventional expression 
of T\ [48]. The solid and the dashed curves in Figs. 1, 
2,6, and 7 are the least-squares fitting results using this 
approximate function. (Ti)h os t is the Cu nuclear spin- 
lattice relaxation time due to the host Cu electron spin 
fluctuation via a hyperfine coupling, n is the impurity- 
induced nuclear spin-lattice relaxation time, i.e., a spatial 
average of l/7i(r) = A(r) 2 Si mp (r, u n ) (r is a distance be- 
tween a nuclear and an impurity-induced moment, A(r) 
is a coupling constant between the nuclear spin and the 
impurity-induced moment, Si mp (r,u) n ) is the impurity- 
induced spin-spin correlation function [9-11], and w n /27r 
is the nuclear resonance frequency) [44,48]. For an iso- 
lated local moment on the impurity site, one obtains 
1/Xi(r) oc A(r) 2 jT (a longitudinal direct dipole coupling 
A(r) 2 cx 1/r 6 , or a two-dimensional RKKY interaction 
A(r) 2 cx 1/r 4 ) [44]. For an impurity-induced moment 
with a local coupling A(r) « A(0), if Si mp (r,ui n ) oc r~ d 
or e _r /£, one obtains p(t) « exp[— (i/ri)™] or exp(-ln 2 i) 
(£ is a correlation length, and a oc d^ 1 is a constant) [45]. 

Instead of including all the direct or indirect nuclear- 
electron interactions, we assume a minimal model with 
a single impurity-induced relaxation and the host Cu 
electron spin relaxation. The recovery curve is ex- 
pressed by p(t) = p(0) exp^^Ti),^] nj(l - c) + 
cexp(— 3£/Ti(rj))] (i is the Cu site index) as a function of 
a lot of time constants of (Ti) host , Ti(ri), Ti(r 2 ), Ti(r 3 ), 
• • •. For c <Cl, the product J7J' ' '] i s approximated by 
« exp[— c / cir(l — exp(— 3t/Ti(r)))], and then the spatial 
integral leads to the stretched exponential function with 



Optimally-doped YBa 2 (Cu I x Zn x ) 3 7 s 
T>69K 



T c =46 K x=0.033 

Satellite 
Main 




100 200 
T(K) 

FIG. 8. Zn-doping effect on the main and the satellite 
63 Cu(2) nuclear spin-lattice relaxation rates 63 (1/tiT) and 
63 (l/TiT)host as functions of temperature for Y123 (a) and 
for Y124 (b). The dashed lines indicate the respective T c 's. T 3 
is the pseudo spin-gap temperature defined by the maximum 
temperature of 63 (l/TiT) host . 



the single time constant T\ [44]. Eventually, only two 
fitting parameters on the time constants, i.e. {T\)h ost 
and Ti are obtained. The distribution of T\ is taken into 
consideration through T\(r), although relations among 
Ti(r), the T 2 distribution, and the broad Cu NQR fre- 
quency spectrum are not straightforward. 

In the actual fitting, we noticed that it is not easy to 
estimate the precise value of 1/^T\)host i n the case 01 
1/ti 3> l/(7i)/ lost . However, Zn-substitution effect on 
each relaxation time is fairly obvious. The slow relax- 
ation of p(t > 1.0 ms) due to heavy doping of Zn as in 
Fig. 6(b) is the experimental fact. Within the present 
model, since p(t) — > p(0) exp[—3t/(Ti)h oat ] at t » n/3, 
then one estimates to decrease l/(Ti)^ ost with Zn doping 
of £=0.033. This estimation is not due to delicate fitting. 
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3. Zn- substitution effect on 63 ti and 63 (Ti) host 

Figure 8 shows the estimated 63 (1/tiT) and 
e3 (l/TiT)host of the satellite and of the main signals 
as functions of temperature for Zn-substituted Y123 (a) 
and Y124 (b). First, the satellite 63 (l/nT) for Y124 at 
T=4.2 K is an order of magnitude higher than that for 
Y123. Second, the temperature dependence of 63 (l/riT) 
of the satellite signal is different between Y123 and Y124. 
For Y124 at lower temperatures below T c , 63 (1/tiT) in- 
creases as temperature is decreased, in contrast to T 
independent or slight decrease of 63 (1/tiT) for Y123. 
Third, the difference in the T dependence of 63 (1/tiT) 
between the main and the satellite signals is marked be- 
low T c for Y123 but above T c for Y124. Fourth, the main 
63 (l/TiT)host above about 100 K decreases with Zn dop- 
ing for both systems. Fifth, the satellite 63 (l/TiT) ftost 
was estimated to be <C 1.0 s _1 K _1 for Y124 with Zn. 
From these results, one can see that the local magnetic 
correlation around Zn is different between undcrdoped 
and optimally doped samples. 

The Cu(2) nuclear spin-lattice relaxation time probes 
the in-plane antiferromagnetic dynamical spin suscepti- 
bility at an NMR/NQR frequency [49-52] and homo- 
geneity of the CuC>2 plane in real space [44,48]. Relation 
between our results on the Cu NQR relaxation and the 
previous 89 Y NMR results [3] is not straightforward. 

6S (1/tiT): The obvious difference in 63 (1/tiT) 
between the main and the satellite signals indi- 
cates an evidence for the Zn-induccd virtual bound 
state [53-55] through the locally enhanced magnetic cor- 
relation [10,11,56,57]. For underdoped Y124, the differ- 
ence in 63 (l/nT) may indicate that Zn-induced "local- 
ized moments" persist both above and below T c [58]. 
For optimally doped Y123, the below- T c difference in 
63 (1/tiT) may indicate that the Zn-induced, generalized 
RKKY oscillation is extended around Zn below T c but 
localized above T c . This is consistent with a prediction 
in the perturbative mode- mode coupling theory [11,56], 
where the spatially extended potential scattering by Zn 
plays a key role [59,60]. The separated frequency of the 
satellite Cu NQR spectrum may also support such a fi- 
nite range effect of Zn impurity potential. According to 
Ref. [11], the 1st nn potential term | V2 j will be larger in 
Y124 than Y123, because the satellite 63 (1/tiT) is larger 
than the main one in Y124 above and below T c . The lo- 
cal electronic state induced by Zn reflects the matrix an- 
tiferromagnetic correlation [10,11,56,57]. The optimally 
doped Y123 still has the underlying antiferromagnetic 
correlation [30,61]. 

63 {1/T! T) host : For Y124, the absence or a small value 
of the satellite 63 (l/Ti7% ost sug gests that the conduc- 
tion electrons in the underdoped system are expelled 
from the vicinity of Zn atoms, or that the pseudo spin-gap 
effect is enhanced near Zn. A remnant of the pseudo spin- 



gap effect can be seen as a decrease of 63 (1/tiT) below 
T=100 K with x=0.005 and 0.01 in Fig. 8(b). For Y123, 
nearly the same 63 (l/TiT)/j OSt between the satellite and 
the main signals indicates that the conduction electrons 
uniformly distribute near and away from Zn both above 
and below T c , in contrast to the underdoped Y124. The 
pure Y124 is a stoichiometric, homogeneous, undcrdoped 
system. Nevertheless, it has some inherent electronic in- 
stability, e.g. magnetic instability, which may lead to 
an inhomogeneous magnetic response. The undcrdoped 
electronic state of Y124 is located more close to the mag- 
netic instability than the optimally doped one of Y123. 

The above- T c decrease in 63 (l/TiT)f lost with Zn doping 
for both systems indicates that Zn suppresses the host an- 
tiferromagnetic spin susceptibility [46]. This is recently 
reproduced by a numerical quantum Monte Carlo sim- 
ulation with the dynamical cluster approximation based 
on the two-dimensional t — J model [62] . 

IV. CONCLUSION 

We succeeded in measuring the Zn-induced satellite 
Cu(2) NQR nuclear spin-lattice relaxation curves in the 
superconductors Y123 and Y124 above and below T c in 
a zero external magnetic field. With doping a small 
amount of Zn, the electron spin correlation near Zn is 
enhanced in the optimally doped system as well as in 
the underdoped one. The optimally doped Y123 still has 
an underlying antiferromagnetic correlation. The local 
electronic state induced by Zn reflects it. 
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